THEORETICAL STUDIES ON THE METHODS OF RECONSTRUCTING PHYLOGENETIC TREES FROM DNA SEQUENCE DATA (MOLECULAR EVOLUTION, HOMINOID EVOLUTION, COMPUTER SIMULATION)

(1) A mathematical theory for computing the probabilities of various nucleotide configurations is developed, and the probability of obtaining the correct phylogenetic tree (model tree) from sequence data is evaluated for six phylogenetic tree-making methods (UPGMA, distance Wagner method, transformed distance method, Fitch-Margoliash's method, maximum parsimony method, and compatibility method). The number of nucleotides (m*) necessary to obtain the correct tree with a probability of 95% is estimated with special reference to the human, chimpanzee, and gorilla divergence. m* is at least 4,200, but the availability of outgroup species greatly reduces m* for all methods except UPGMA. m* increases if transitions occur more frequently than transversions as in the case of mitochondrial DNA. (2) A new tree-making method called the neighbor-joining method is proposed. This method is applicable either for distance data or character state data. Computer simulation has shown that the neighbor-joining method is generally better than UPGMA, Farris' method, Li's method, and modified Farris method on recovering the true topology when distance data are used. A related method, the simultaneous partitioning method, is also discussed. (3) The maximum likelihood (ML) method for phylogeny reconstruction under the assumption of both constant and varying evolutionary rates is studied, and a new algorithm for obtaining the ML tree is presented. This method gives a tree similar to that obtained by UPGMA when constant evolutionary rate is assumed, whereas it gives a tree similar to that obtained by the maximum parsimony tree and the neighbor-joining method when varying evolutionary rate is assumed. ^

Detailed analysis of sequence changes occurring during vlsE antigenic variation in the mouse model of Borrelia burgdorferi infection.

Lyme disease Borrelia can infect humans and animals for months to years, despite the presence of an active host immune response. The vls antigenic variation system, which expresses the surface-exposed lipoprotein VlsE, plays a major role in B. burgdorferi immune evasion. Gene conversion between vls silent cassettes and the vlsE expression site occurs at high frequency during mammalian infection, resulting in sequence variation in the VlsE product. In this study, we examined vlsE sequence variation in B. burgdorferi B31 during mouse infection by analyzing 1,399 clones isolated from bladder, heart, joint, ear, and skin tissues of mice infected for 4 to 365 days. The median number of codon changes increased progressively in C3H/HeN mice from 4 to 28 days post infection, and no clones retained the parental vlsE sequence at 28 days. In contrast, the decrease in the number of clones with the parental vlsE sequence and the increase in the number of sequence changes occurred more gradually in severe combined immunodeficiency (SCID) mice. Clones containing a stop codon were isolated, indicating that continuous expression of full-length VlsE is not required for survival in vivo; also, these clones continued to undergo vlsE recombination. Analysis of clones with apparent single recombination events indicated that recombinations into vlsE are nonselective with regard to the silent cassette utilized, as well as the length and location of the recombination event. Sequence changes as small as one base pair were common. Fifteen percent of recovered vlsE variants contained "template-independent" sequence changes, which clustered in the variable regions of vlsE. We hypothesize that the increased frequency and complexity of vlsE sequence changes observed in clones recovered from immunocompetent mice (as compared with SCID mice) is due to rapid clearance of relatively invariant clones by variable region-specific anti-VlsE antibody responses.

Molecular basis for class V beta-tubulin effects on microtubule assembly and paclitaxel resistance.

Vertebrates produce at least seven distinct beta-tubulin isotypes that coassemble into all cellular microtubules. The functional differences among these tubulin isoforms are largely unknown, but recent studies indicate that tubulin composition can affect microtubule properties and cellular microtubule-dependent behavior. One of the isotypes whose incorporation causes the largest change in microtubule assembly is beta5-tubulin. Overexpression of this isotype can almost completely destroy the microtubule network, yet it appears to be required in smaller amounts for normal mitotic progression. Moderate levels of overexpression can also confer paclitaxel resistance. Experiments using chimeric constructs and site-directed mutagenesis now indicate that the hypervariable C-terminal region of beta5 plays no role in these phenotypes. Instead, we demonstrate that two residues found in beta5 (Ser-239 and Ser-365) are each sufficient to inhibit microtubule assembly and confer paclitaxel resistance when introduced into beta1-tubulin; yet the single mutation of residue Ser-239 in beta5 eliminates its ability to confer these phenotypes. Despite the high degree of conservation among beta-tubulin isotypes, mutations affecting residue 365 demonstrate that amino acid substitutions can be context sensitive; i.e. an amino acid change in one isotype will not necessarily produce the same phenotype when introduced into a different isotype. Modeling studies indicate that residue Cys-239 of beta1-tubulin is close to a highly conserved Cys-354 residue suggesting the possibility that disulfide formation could play a significant role in the stability of microtubules formed with beta1- but not with beta5-tubulin.

Mutations in the cofilin partner Aip1/Wdr1 cause autoinflammatory disease and macrothrombocytopenia.

A pivotal mediator of actin dynamics is the protein cofilin, which promotes filament severing and depolymerization, facilitating the breakdown of existing filaments, and the enhancement of filament growth from newly created barbed ends. It does so in concert with actin interacting protein 1 (Aip1), which serves to accelerate cofilin's activity. While progress has been made in understanding its biochemical functions, the physiologic processes the cofilin/Aip1 complex regulates, particularly in higher organisms, are yet to be determined. We have generated an allelic series for WD40 repeat protein 1 (Wdr1), the mammalian homolog of Aip1, and report that reductions in Wdr1 function produce a dramatic phenotype gradient. While severe loss of function at the Wdr1 locus causes embryonic lethality, macrothrombocytopenia and autoinflammatory disease develop in mice carrying hypomorphic alleles. Macrothrombocytopenia is the result of megakaryocyte maturation defects, which lead to a failure of normal platelet shedding. Autoinflammatory disease, which is bone marrow-derived yet nonlymphoid in origin, is characterized by a massive infiltration of neutrophils into inflammatory lesions. Cytoskeletal responses are impaired in Wdr1 mutant neutrophils. These studies establish an essential requirement for Wdr1 in megakaryocytes and neutrophils, indicating that cofilin-mediated actin dynamics are critically important to the development and function of both cell types.

Eomesodermin, a target gene of Pou4f2, is required for retinal ganglion cell and optic nerve development in the mouse.

The mechanisms regulating retinal ganglion cell (RGC) development are crucial for retinogenesis and for the establishment of normal vision. However, these mechanisms are only vaguely understood. RGCs are the first neuronal lineage to segregate from pluripotent progenitors in the developing retina. As output neurons, RGCs display developmental features very distinct from those of the other retinal cell types. To better understand RGC development, we have previously constructed a gene regulatory network featuring a hierarchical cascade of transcription factors that ultimately controls the expression of downstream effector genes. This has revealed the existence of a Pou domain transcription factor, Pou4f2, that occupies a key node in the RGC gene regulatory network and that is essential for RGC differentiation. However, little is known about the genes that connect upstream regulatory genes, such as Pou4f2 with downstream effector genes responsible for RGC differentiation. The purpose of this study was to characterize the retinal function of eomesodermin (Eomes), a T-box transcription factor with previously unsuspected roles in retinogenesis. We show that Eomes is expressed in developing RGCs and is a mediator of Pou4f2 function. Pou4f2 directly regulates Eomes expression through a cis-regulatory element within a conserved retinal enhancer. Deleting Eomes in the developing retina causes defects reminiscent of those in Pou4f2(-/-) retinas. Moreover, myelin ensheathment in the optic nerves of Eomes(-/-) embryos is severely impaired, suggesting that Eomes regulates this process. We conclude that Eomes is a crucial regulator positioned immediately downstream of Pou4f2 and is required for RGC differentiation and optic nerve development.

Sp1 up-regulates cAMP-response-element-binding protein expression during retinoic acid-induced mucous differentiation of normal human bronchial epithelial cells.

CREB [CRE (cAMP-response element)-binding protein] is an important transcription factor that is differentially regulated in cells of various types. We recently reported that RA (retinoic acid) rapidly activates CREB without using RARs (RA receptors) or RXRs (retinoid X receptors) in NHTBE cells (normal human tracheobronchial epithelial cells). However, little is known about the role of RA in the physiological regulation of CREB expression in the early mucous differentiation of NHTBE cells. In the present study, we report that RA up-regulates CREB gene expression and that, using 5'-serial deletion promoter analysis and mutagenesis analyses, two Sp1 (specificity protein 1)-binding sites located at nt -217 and -150, which flank the transcription initiation site, are essential for RA induction of CREB gene transcription. Furthermore, we found that CREs located at nt -119 and -98 contributed to basal promoter activity. Interestingly, RA also up-regulated Sp1 in a time- and dose-dependent manner. Knockdown of endogenous Sp1 using siRNA (small interfering RNA) decreased RA-induced CREB gene expression. However, the converse was not true: knockdown of CREB using CREB siRNA did not affect RA-induced Sp1 gene expression. We conclude that RA up-regulates CREB gene expression during the early stage of NHTBE cell differentiation and that RA-inducible Sp1 plays a major role in up-regulating human CREB gene expression. This result implies that co-operation of these two transcription factors plays a crucial role in mediating early events of normal mucous cell differentiation of bronchial epithelial cells.

Rtr1 is the Saccharomyces cerevisiae homolog of a novel family of RNA polymerase II-binding proteins.

Cells must rapidly sense and respond to a wide variety of potentially cytotoxic external stressors to survive in a constantly changing environment. In a search for novel genes required for stress tolerance in Saccharomyces cerevisiae, we identified the uncharacterized open reading frame YER139C as a gene required for growth at 37 degrees C in the presence of the heat shock mimetic formamide. YER139C encodes the closest yeast homolog of the human RPAP2 protein, recently identified as a novel RNA polymerase II (RNAPII)-associated factor. Multiple lines of evidence support a role for this gene family in transcription, prompting us to rename YER139C RTR1 (regulator of transcription). The core RNAPII subunits RPB5, RPB7, and RPB9 were isolated as potent high-copy-number suppressors of the rtr1Delta temperature-sensitive growth phenotype, and deletion of the nonessential subunits RPB4 and RPB9 hypersensitized cells to RTR1 overexpression. Disruption of RTR1 resulted in mycophenolic acid sensitivity and synthetic genetic interactions with a number of genes involved in multiple phases of transcription. Consistently, rtr1Delta cells are defective in inducible transcription from the GAL1 promoter. Rtr1 constitutively shuttles between the cytoplasm and nucleus, where it physically associates with an active RNAPII transcriptional complex. Taken together, our data reveal a role for members of the RTR1/RPAP2 family as regulators of core RNAPII function.

A loss-of-function variant of PTPN22 is associated with reduced risk of systemic lupus erythematosus.

A gain-of-function R620W polymorphism in the PTPN22 gene, encoding the lymphoid tyrosine phosphatase LYP, has recently emerged as an important risk factor for human autoimmunity. Here we report that another missense substitution (R263Q) within the catalytic domain of LYP leads to reduced phosphatase activity. High-resolution structural analysis revealed the molecular basis for this loss of function. Furthermore, the Q263 variant conferred protection against human systemic lupus erythematosus, reinforcing the proposal that inhibition of LYP activity could be beneficial in human autoimmunity.

Identification and characterization of the herpes simplex virus type 2 gene encoding the essential capsid protein ICP32/VP19c.

We describe the characterization of the herpes simplex virus type 2 (HSV-2) gene encoding infected cell protein 32 (ICP32) and virion protein 19c (VP19c). We also demonstrate that the HSV-1 UL38/ORF.553 open reading frame (ORF), which has been shown to specify a viral protein essential for capsid formation (B. Pertuiset, M. Boccara, J. Cebrian, N. Berthelot, S. Chousterman, F. Puvian-Dutilleul, J. Sisman, and P. Sheldrick, J. Virol. 63: 2169-2179, 1989), must encode the cognate HSV type 1 (HSV-1) ICP32/VP19c protein. The region of the HSV-2 genome deduced to contain the gene specifying ICP32/VP19c was isolated and subcloned, and the nucleotide sequence of 2,158 base pairs of HSV-2 DNA mapping immediately upstream of the gene encoding the large subunit of the viral ribonucleotide reductase was determined. This region of the HSV-2 genome contains a large ORF capable of encoding two related 50,538- and 49,472-molecular-weight polypeptides. Direct evidence that this ORF encodes HSV-2 ICP32/VP19c was provided by immunoblotting experiments that utilized antisera directed against synthetic oligopeptides corresponding to internal portions of the predicted polypeptides encoded by the HSV-2 ORF or antisera directed against a TrpE/HSV-2 ORF fusion protein. The type-common immunoreactivity of the two antisera and comparison of the primary amino acid sequences of the predicted products of the HSV-2 ORF and the equivalent genomic region of HSV-1 provided evidence that the HSV-1 UL38 ORF encodes the HSV-1 ICP32/VP19c. Analysis of the expression of the HSV-1 and HSV-2 ICP32/VP19c cognate proteins indicated that there may be differences in their modes of synthesis. Comparison of the predicted structure of the HSV-2 ICP32/VP19c protein with the structures of related proteins encoded by other herpes viruses suggested that the internal capsid architecture of the herpes family of viruses varies substantially.

A Cbfa1-dependent genetic pathway controls bone formation beyond embryonic development.

The molecular mechanisms controlling bone extracellular matrix (ECM) deposition by differentiated osteoblasts in postnatal life, called hereafter bone formation, are unknown. This contrasts with the growing knowledge about the genetic control of osteoblast differentiation during embryonic development. Cbfa1, a transcriptional activator of osteoblast differentiation during embryonic development, is also expressed in differentiated osteoblasts postnatally. The perinatal lethality occurring in Cbfa1-deficient mice has prevented so far the study of its function after birth. To determine if Cbfa1 plays a role during bone formation we generated transgenic mice overexpressing Cbfa1 DNA-binding domain (DeltaCbfa1) in differentiated osteoblasts only postnatally. DeltaCbfa1 has a higher affinity for DNA than Cbfa1 itself, has no transcriptional activity on its own, and can act in a dominant-negative manner in DNA cotransfection assays. DeltaCbfa1-expressing mice have a normal skeleton at birth but develop an osteopenic phenotype thereafter. Dynamic histomorphometric studies show that this phenotype is caused by a major decrease in the bone formation rate in the face of a normal number of osteoblasts thus indicating that once osteoblasts are differentiated Cbfa1 regulates their function. Molecular analyses reveal that the expression of the genes expressed in osteoblasts and encoding bone ECM proteins is nearly abolished in transgenic mice, and ex vivo assays demonstrated that DeltaCbfa1-expressing osteoblasts were less active than wild-type osteoblasts. We also show that Cbfa1 regulates positively the activity of its own promoter, which has the highest affinity Cbfa1-binding sites characterized. This study demonstrates that beyond its differentiation function Cbfa1 is the first transcriptional activator of bone formation identified to date and illustrates that developmentally important genes control physiological processes postnatally.

Construction of improved temperature-sensitive and mobilizable vectors and their use for constructing mutations in the adhesin-encoding acm gene of poorly transformable clinical Enterococcus faecium strains.

Inactivation by allelic exchange in clinical isolates of the emerging nosocomial pathogen Enterococcus faecium has been hindered by lack of efficient tools, and, in this study, transformation of clinical isolates was found to be particularly problematic. For this reason, a vector for allelic replacement (pTEX5500ts) was constructed that includes (i) the pWV01-based gram-positive repAts replication region, which is known to confer a high degree of temperature intolerance, (ii) Escherichia coli oriR from pUC18, (iii) two extended multiple-cloning sites located upstream and downstream of one of the marker genes for efficient cloning of flanking regions for double-crossover mutagenesis, (iv) transcriptional terminator sites to terminate undesired readthrough, and (v) a synthetic extended promoter region containing the cat gene for allelic exchange and a high-level gentamicin resistance gene, aph(2'')-Id, to distinguish double-crossover recombination, both of which are functional in gram-positive and gram-negative backgrounds. To demonstrate the functionality of this vector, the vector was used to construct an acm (encoding an adhesin to collagen from E. faecium) deletion mutant of a poorly transformable multidrug-resistant E. faecium endocarditis isolate, TX0082. The acm-deleted strain, TX6051 (TX0082Deltaacm), was shown to lack Acm on its surface, which resulted in the abolishment of the collagen adherence phenotype observed in TX0082. A mobilizable derivative (pTEX5501ts) that contains oriT of Tn916 to facilitate conjugative transfer from the transformable E. faecalis strain JH2Sm::Tn916 to E. faecium was also constructed. Using this vector, the acm gene of a nonelectroporable E. faecium wound isolate was successfully interrupted. Thus, pTEX5500ts and its mobilizable derivative demonstrated their roles as important tools by helping to create the first reported allelic replacement in E. faecium; the constructed this acm deletion mutant will be useful for assessing the role of acm in E. faecium pathogenesis using animal models.

Translational regulation of nuclear gene COX4 expression by mitochondrial content of phosphatidylglycerol and cardiolipin in Saccharomyces cerevisiae.

Previous results indicated that translation of four mitochondrion-encoded genes and one nucleus-encoded gene (COX4) is repressed in mutants (pgs1Delta) of Saccharomyces cerevisiae lacking phosphatidylglycerol and cardiolipin. COX4 translation was studied here using a mitochondrially targeted green fluorescence protein (mtGFP) fused to the COX4 promoter and its 5' and 3' untranslated regions (UTRs). Lack of mtGFP expression independent of carbon source and strain background was established to be at the translational level. The translational defect was not due to deficiency of mitochondrial respiratory function but was rather caused directly by the lack of phosphatidylglycerol and cardiolipin in mitochondrial membranes. Reintroduction of a functional PGS1 gene under control of the ADH1 promoter restored phosphatidylglycerol synthesis and expression of mtGFP. Deletion analysis of the 5' UTR(COX4) revealed the presence of a 50-nucleotide fragment with two stem-loops as a cis-element inhibiting COX4 translation. Binding of a protein factor(s) specifically to this sequence was observed with cytoplasm from pgs1Delta but not PGS1 cells. Using HIS3 and lacZ as reporters, extragenic spontaneous recessive mutations that allowed expression of His3p and beta-galactosidase were isolated, which appeared to be loss-of-function mutations, suggesting that the genes mutated may encode the trans factors that bind to the cis element in pgs1Delta cells.

Functional analysis of the amine substrate specificity domain of pepper tyramine and serotonin N-hydroxycinnamoyltransferases.

Pepper (Capsicum annuum) serotonin N-hydroxycinnamoyltransferase (SHT) catalyzes the synthesis of N-hydroxycinnamic acid amides of serotonin, including feruloylserotonin and p-coumaroylserotonin. To elucidate the domain or the key amino acid that determines the amine substrate specificity, we isolated a tyramine N-hydroxycinnamoyltransferase (THT) gene from pepper. Purified recombinant THT protein catalyzed the synthesis of N-hydroxycinnamic acid amides of tyramine, including feruloyltyramine and p-coumaroyltyramine, but did not accept serotonin as a substrate. Both the SHT and THT mRNAs were found to be expressed constitutively in all pepper organs. Pepper SHT and THT, which have primary sequences that are 78% identical, were used as models to investigate the structural determinants responsible for their distinct substrate specificities and other enzymatic properties. A series of chimeric genes was constructed by reciprocal exchange of DNA segments between the SHT and THT cDNAs. Functional characterization of the recombinant chimeric proteins revealed that the amino acid residues 129 to 165 of SHT and the corresponding residues 125 to 160 in THT are critical structural determinants for amine substrate specificity. Several amino acids are strongly implicated in the determination of amine substrate specificity, in which glycine-158 is involved in catalysis and amine substrate binding and tyrosine-149 plays a pivotal role in controlling amine substrate specificity between serotonin and tyramine in SHT. Furthermore, the indisputable role of tyrosine is corroborated by the THT-F145Y mutant that uses serotonin as the acyl acceptor. The results from the chimeras and the kinetic measurements will direct the creation of additional novel N-hydroxycinnamoyltransferases from the various N-hydroxycinnamoyltransferases found in nature.

Evolution of sensory complexity recorded in a myxobacterial genome.

Myxobacteria are single-celled, but social, eubacterial predators. Upon starvation they build multicellular fruiting bodies using a developmental program that progressively changes the pattern of cell movement and the repertoire of genes expressed. Development terminates with spore differentiation and is coordinated by both diffusible and cell-bound signals. The growth and development of Myxococcus xanthus is regulated by the integration of multiple signals from outside the cells with physiological signals from within. A collection of M. xanthus cells behaves, in many respects, like a multicellular organism. For these reasons M. xanthus offers unparalleled access to a regulatory network that controls development and that organizes cell movement on surfaces. The genome of M. xanthus is large (9.14 Mb), considerably larger than the other sequenced delta-proteobacteria. We suggest that gene duplication and divergence were major contributors to genomic expansion from its progenitor. More than 1,500 duplications specific to the myxobacterial lineage were identified, representing >15% of the total genes. Genes were not duplicated at random; rather, genes for cell-cell signaling, small molecule sensing, and integrative transcription control were amplified selectively. Families of genes encoding the production of secondary metabolites are overrepresented in the genome but may have been received by horizontal gene transfer and are likely to be important for predation.

Insights into transcription enhancer factor 1 (TEF-1) activity from the solution structure of the TEA domain.

Transcription enhancer factor 1 is essential for cardiac, skeletal, and smooth muscle development and uses its N-terminal TEA domain (TEAD) to bind M-CAT elements. Here, we present the first structure of TEAD and show that it is a three-helix bundle with a homeodomain fold. Structural data reveal how TEAD binds DNA. Using structure-function correlations, we find that the L1 loop is essential for cooperative loading of TEAD molecules on to tandemly duplicated M-CAT sites. Furthermore, using a microarray chip-based assay, we establish that known binding sites of the full-length protein are only a subset of DNA elements recognized by TEAD. Our results provide a model for understanding the regulation of genome-wide gene expression during development by TEA/ATTS family of transcription factors.